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I. UCTION 


Often, in the design and production of complex system such as ccmputer 
network, good system reliability is a critical requirement. Network reliability is 
rarely treated explicitly and quantitatively, however, largely because of problems 
of computational complexity and because few efficient computer programs exist 
for computing network reliability. Only in the last few years have research and 


development in this area provided usable results. 


We wish, as efficiently as possible, to determine the network reliability for 
specific class of networks. The purpose of this thesis is to develop new reduction 
and decomposition techniques for computing network reliability and to show how 


computational complexity can be decreased by using these techniques. 


The network model which is used in this thesis may be thousht of as a 
communication network with directed edges, allowing only one way communication 
from tail vertex to head vertex. Additionally, we impose the restriction that the 
network have no directed cycles. Such a model may be appropriate in hierarchical 
communication networks such as command networks. Also, such networks may 


arise as subproblems while computing the reliability of general directed networks, 


Agrawal [1982]. 


Networks can be divided into various classes, e.g., directed, undirected, cyclic, 
acyclic, etc. A large amount of literature exists on each class. (See Agrawal 
[1982], Ball and Provan [1981], Satyanarayana [1982], Hagstrom [1980] and 
Satyanarayana and Prabahkar [1978] for directed networks; see Ball and Provan 
[1981], Johnson [1982], Satyanarayana [1982], Satyanarayana and Chang [1981], 


Satyanaryana and Wood (1982], Valiant ¡1978], and Wood [1982] for undirected 





networks; see Agrawal [1982] and Satyanarayara and Prabhakar [1978] for cyclic 
versus acyclic directed networks.) The specific model we will consider is the 
directed acyclic network, denoted D(V,E) where V is the set of vertices and £ is 
the set of edges. Each edge is represented by an ordered pair of vertices e =(u,v). 
The ordering implies that the flow of communication can move along e from u to v 
only. 

Each edge e, in the network functions with probability p, and fails with 
probability بو‎ = 1—p,. All edge failures are assumed to be independent events. For 
simplicity in the initial discussion, vertices are assumed to be completely reliable. 


However, the extension to unreliable vertices will be handled briefly. 


Vertex u can send communication to vertex v in D if there is a directed path 
of working edges from u tov. One of the vertices s is designated as the "source" 
or the "root". With respect to the source, there are three different measures of 
reliability which are usually studied in directed networks: 

(1) The source-to-terminal (ST) reliability, A,,(D), is the probability that s can 
send communication to a specified vertex f in J. 

(2) The source-to-all-terminal (SAT) reliability, R,y(D), is the probability that s can 
send communication to all vertices of D. 

(3) The source-to-K-terminal (SKT) reliability, R,g(D), is the probability that s can 
send communication to every vertex in a specified set K (these vertices are 
referred to as "K-vertices") with, by convention, s eK. 

Of course, (1) and (2) are special cases of (3). We consider the general case (2) in 

this tnesis. 

Generally, the of problem computing network reliability is NP-hard. Certain 
reduction and decomposition schemes exist which reduce the size oí the problem 
while preserving network reliability, e.g., Agrawal[1982]. Shogan[:978] and Wood 


[1982], and, for special classes of networks, it is thereby pcssible to compute 





reliability in polynomial time. In general, however, the solution time grows 
exponentially in the size of the network. In this thesis, we develop a method for 
computing reliability which, in general, is exponential in the size of D, but which is 
polynomial when the number of non-K-vertices, i.e., .|Y-X |, is fixed. 

The rest of this thesis is outlined as follows. 

Chapter 2 contains a necessary definitions and graph theoretic concepts. 

Chapter 3 develops reduction and decomposition techniques to reduce the size 
of a network in order to decrease computation time. Some of these techniques are 
new. 

Chapter 4 shows how to compute network reliability exactly when the number 
of non-K-vertices is moderate. A truncated version of this technique is used to 
obtain a lower bound on reliability when the number of non-K-vertices is too large. 
We demonstrate that tnis lower bound will usually be "zood." 

Chapter 5 describes an algorithm to compute R,x(D) or find a lower bound on 
R,g(D), and its implementation. The program is coded in FORTRAN with one main 
routine and 18 subroutines totaling about 1600 lines of code. 


Chapter 6 gives computational experience for the algorithm run on an IBM 


3033AP system /370 with the FORTRAN H (Extended) compiler. 


Chapter 7 is a conclusion and gives suggestions for further research 








In this chapter, we give the definitions and notation used throuzhout this 


thesis, and introduce basic network reliability theory. 


A. DEFINITIONS AND NOTATION 


A directed network D=(V,£) comprises two finite sets: V is the set of vertices 
and E is the set of edges. Each edge eel corresponds to an ordered pair of 
vertices in the directed network. That is, e z(u,v) where u,vcV. The vertex u is 
called the tail vertex of edge e and vertex v is called the head vertex oí edge e. 
For any two vertices u and v in D, if there exists an edge e -(u,v), then u and v 
are said to be adjacent and edge e is incident into v and incident out of of u. The 
indesree of a vertex v, denoted indeg(v), is the number of edges incident into v 
and the outdegree of a vertex v, denoted outdeg(v), is the number of edges 
incident out of v. Two vertices are connected (or communicate") if there exists a 
sequence of vertices and edges of the form vo, (vov), vi, (v1.2)... (Um-1: Um). Ym- 


This sequence is a path of length m. If the vg-v,, the path is a cycle. 


Consider a undirected graph G(V,E) fcrraed by ignoring the ordering of the 
edges in a directed graph D(V,E). G is called the underlying graph of D. D is 
connected if its underlying graph G is connected. G is connected if there is an 
undirected path between all pairs of vertices. The vertex connectivity of D is the 
minimum number of vertices which must be deleted from G (along with adjacent 
edges, of course) in order to disconnect G or to create a trivial graph with only one 


vertex. D is separable if its vertex connectivity is one. Otherwise D is 





nonseparable. In a separable graph, a vertex whose deletion disconnects the graph 


is called a cut-vertex 


Merging a vertex v into another vertex w means removing all edges between v 
and w and forming a single vertex uv’ by coalescing v and w such that any edge 
incident into (out of) v or w is now incident into (out of) v'. Merging v and w is 


sometimes also referred to as "identifying" v with w. 


The SKT reliability of a directed network D, denoted K,&(D), is the probability 
that all K-vertices in D are connected by path of working edges and vertices from 
source vertex s. The probability ihat edge e; works is p; = 1-97;. For now, we 


assume that all vertices work with probability one. 


B. NETWORK THEORY AND SURVEY 
1. NP-Completeness 

We briefly discuss the complexity of the SKT reliability problem in 
terms of the theory of NP-completeness and in terms of practical computation. 
More detail can be found in Garey and Johnson | 1979], Ball and Provan| 1981]. 

Algorithms to solve various problems can be broadly classified into two 
categories, polynomial-time and exponential-time algorithms. An algorithm is a 
polynomial-time algorithm if for a problem of size n, its running time is bounded 
by a polynomial in n. Any algorithm that is not a polynomial-time algorithm is an 
exponential-time algorithm. In combinatorics, a decision problem P is said to 
belong to the class NP if given a tentative soiution, the validity of the solution can 
be checked in polynomial-time. P is "NP-complete" if it is equivalent to the 
satisfiability problem of a conjunctive normal form boolean expression. A problem 


is NP-hard if it is at least as hard as an NP-complete problem. 
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A problem is #P-complete (number-P complete) if its solution is equivalent to 
counting (not enumerating) the number of solutions to an NP-complete problem. 
Any #P-complete problem is NP-hard. 
2. PBrablemCamplesity 
Ball and Provan [1981] show that the ST reliability problem in an acyclic 
directed network is a #P-complete problem. This implies that the SKT problem is 
also #P-complete since problem instances restricted to |K|=2 are #P-complete. 
Even finding an -approximation to SKT reliability is an #P-complete problem. 
| In practice, these complexity results imply that a generel problem of 
large size cannot be solved analytically in a reasonable amount of time. In fact, 
the time taken grows exponentialy with size of the network. As will be seen, 
however, practical problems of fairly large size can be handled with the techniques 
developed here. Furthermore, even larger problems can be handled if one is 
willing to accept a "good" lower bound on exact reliability. 
J. | BRedurticn 
Many reduction schemes have been developed to help solve various 
network reliability problems. These reductions reduce the size of the network in 
polynomial time while preserving reliability, i.e., they create a new, smaller, 
network G' from the original network G (G may directed or undirected and any 
number of definitions of reliability may be used.) such that R(G)=M R(G') where M 
is a known constant factor. In this way, a polynomial amount cf work is expended 
in order to save, it is hoped, an exponential amount of work during later, more 
general computations. Ín some special, but often practical cases, it is possible to 
coinpletely reduce a network and effectively compute its reliability in polynomial 
time. 
The reductions describe below give a flavor for the different types of 
reductions which have been developed. Parallel and series reductions are well- 


known schemes so they will not 5e discussed. The reductions described assume 
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that vertices are completely reliable. However, most can be extended to hancle 
unreliable vertices. 
(1) Polygon-To-Chain Reductions: A chain C in an undirected graph G is defined as 
a alternating sequence of vertices and edges beginning and ending with vertices 
such that all internal vertices have degree 2. Two parallel chains form a polygon 
and this polygon can be replaced by a chain in a reduction that preserves 
reliability in the K-terminal reliability problem, the analog to the SKT reliability 
problem in an undirected network. This reduction and its implications are 
developed by Wood [1982]. 
(2) 2-Neighbor Vertex Reduction: For the SKT reliability problem. any vertex u 
that has exactly two neighbors v and w can be eliminated from the network if 
ug Kor ju,v,wiek anduss. 
(3) Parallel-Chain Reduction: Suppose E, 15 (re set of all edges in a directed 
network D whose endpoints are v; and v;. In D, a chain G is an alternating 
sequence of vertices and non-empty sets of edges where the internal vertices are 
e-neighbor vertices in D. Two chains C, and C; with identical endpoints are 
parallel chains. Parallel chains can be reduced to a single chain in the SKT 
reliability problem. Reductions (2) and (3) are developed by Agrawal [1932]. 
4. Decomposition 

Sometimes it is possible to divide a given graph into a number of 
subgraphs, such that each subgraph can be analyzed separately and the results 
combined to obtain the reliability of the whole graph. We briefly describe some of 
these methods based on vertex connectivity. 
(1) 2-Connected Graph Decomposition: Let G be a 2-connected graph with a pair of 
vertices {u,u} whose deletion disconnects the graph. These vertices are called a 
“separating pair.” This separating pair partitions G into two subgraphs G, and Gg 
such that G,, Ge =G, Gi(\Gz =)u,v}, £,4#¢ and fed. It has been shown by Wood 


[1982] that, in the K-terminal reliability problem, each subgraph defines a chain 
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through a set of reliability computations on the subgraphs. The subgraphs may be 
replaced by the chains so as to preserve reliability. Agrawal [1982] has proven the 
analogous decomposition scheme for the SKT reliability problem on directed 
networks. (2) Strong-connected Decomposition: A strong component of a directed 
graph D is a maximal subgraph of D such that there exists a directed path from 
every vertex to every other vertex in the subgraph. When computing R,x(D), each 
strong component can be considered separately. Computing each componen: 
R,k(Dj) and multipling Rx(D) is whole graph reliability A,x(D). This 
decomposition has been developed by Shogan [1973] and Agrawal [ 1982]. 

Other decomposition schemes exist. Decomposition in 3-connected 
graphs is often possible, Rosenthal [1974]. Decomposition in 1-connected graphs 


will be discussed in chapter 3. 
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111. REDUCTION AND DECOMPOSTTION 


Any algorithm to solve a general network reliability problem has a running 
time which is exponential in the size of problem. If by some method, the size of the 
problem can be reduced, immense computational savings can be made. Therefore, 
in this chapter, we introduce reduction and decomposition schemes which wil! 
often reduce the size of the problem while always preserving reliability. Some of 
these techniques are specific to the SKT reliability problem in acyclic directed 


networks while others are of more general applicability. 


Reduction sckemes reduce a parameter of the size of the problem, usually 
edges or vertices, in order to decrease computational complexity while preserving 
network reliability. On the other hand, decomposition schemes split the whole 
graph into subgraphs, computes subgraphs reliabilities separately and then 


combines these reliabilities in some way to compute overall network reliability. 


A REDUCTIONS 


In this section, we discuss various reduction schemes which reduce the size of 
D and, consequently, reduce the complexity of computing fsx(D). Under any 
reduction, reliability of the network remains invariant up to a known, constant 
factor, i.e., R,x(D) = M'R,œ(D') where H is the known factor, K is the new set of 
K-vertices, and D' is the reduced network. In the following illustrations, all vertices 


ere represented as circles with K-vertices being darkened or shaded. 


The parallel and series reductions are standard, i.e., weil-known and widely 


used. The three reduction schemes, 
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"indegree-one K-vertex contraction", "non-K-vertex deletion" and "neck vertex to 
K-vertex reduction" are newly developed reductions. 
1. Parallel Edges Reduction 
In a directed network, two parallel edges e,=(u,v) and e, =(u,v) may 
be replaced a single edge e, z(u,v) while preserving reliability if p. < و2‎ +۵ —PaDbo 
and M=1. 
2. Series Edges Reduction 
In a directed network, a pair of edges eg=(u,v) and e, =(u,w) such that 
indeg(v )=outdeg(u)=1 are series edges. If ue K, edges eg and ey can be replaced 
by single edge e, =(u,w) while preserving reliability ifp, =p,p, and W=1. 
3. e- ver 
Let e; z(u,v) be en edge in D such that veX and indeg(v)=1. In the SKT 
reliability prcblem, communication from s to vertex veK can only only occur via 
edge e, so-this edge must work if the network is to ar Consequently, edge e; 


can be contracted such that F4 (D) 2 M- F,(D') and Mop 
where D'=| V-u -v +w, E-e? 
w = ulju 


len K-vutw,itusg K 
^| K-v-u-ctw,ifuecK. 


The validity of this reduction is easily showa. Let e;(u,v) be the edge 
described above in the acyclic network D. Let Ff; denote the event that e; is 
working, and let F; denote the complementary event. Since R,x(D) is just a 


probability, the rules of conditional probability can be applied to obtain 


Rax(D) = Pi Rsx(D) Fi)+9: Rsx(D | Fi) 
= p¡Pex(D| Fi since Ayx(D|F,)=0 
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= Pi Rsx(D') 


since no new paths are created by contraction. 
4. -K- 8 


In the SKT reliability problem, an outdegree zero non-K-vertex can be 
deleted from the network without changing network reliability since such a vertex 
is irrelevant. In this case Ryx(D) = MH Ryx(D—wv) and M=1. 

9. zu 

Suppose u and v are two vertices in D such that ueK, vs K and all 
paths from s to u include v. v is called a neck vertex. It may or may not be a 
cut-vertex. Since s must communicate with v before it communicates with u, 
Rex(D) = M Rs(x+v)(D) where M=1. Thus, the neck vertex to K-vertex reduction 
simply changes K to K+v. See Figures 3.1a and 3.ib. (All figures are given in 
Appendix A.) 

The neck vertex to K-vertex is termed a "reduction" because it reduces 
the parameter |V-K| in D. As will be seen in chapter 4, the general algorithm for 
computing A,k(D) runs in time proportional to 2!V-£! so this is a reduction 
according to our definition. Furthermore, the neck vertex to K-vertex reduction 
will often allow other reductions or decomposition to take place. 

6. Other Reductions 

Numercus other reductions exist as exemplified in chapter 2, section 

B3. However, none of these reductions has been implemented because of their 


complexity, and therefore, no other reductions will be described here. 


B. DECOMPOSITION 


Sometimes it is possible to decompose a given network J into a number of 
subnetworks. Each subnetwork can then be analyzed separately and the results 


combined in some way to obtain overall network reliability. 
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Decomposition can result in substantial savings because, in general, the time 
to analyze the whole system is exponential in size of the system. Effort needed to 
compute A,x(D) using decomposition is the sum of efforts needed for each 
subproblem plus the effort, normally polynomial in the size of D, in performing the 
decomposition and recombining the result. Without decomposition, the necessary 
work is rough!y the product of the effort needed to solve the individual problems. 

1. Cut-Vertex Decompnasition 

A connected, undirected, network G(V,E) has a cut-vertex v if there 
exist components (subgraphs) C, and C2 such that C,;(\C2=u and C,;/Ce= G. Any 
communication from one compcenent to the other must pass through v. In this 
case we aiso say that v separates C, and Cs. A network which has a cut-vertex is 
"separable" and one which has none is called "norseparable.” A directed network D 
is separable if its underlying graph G is separable. 

Let D have nonseparable components C,,C5,..., C4. Any component C; 
which has only one cut-vertex is called a pendanf component. If a pendant 
component, G, contains no K-vertices except possibly its cut-vertex, then it is 
completely irrelevant with respect to reliability since none of its edges are in any 
path from the source to a K-vertex. Thus, this component, excluding its cut-vertex, 
can be deleted from the network while preserving network reliability in the SKT 
reliability proolem. This type of deletion may be recursively applied urtil only 
relevant components remain. 

Cut-vertices are not efficiently found in a directed network if each edge 
is represented in one direction, e.g., if all edges incident into each vertex are 
stored contiguously. However, if edges are stored both by tail vertex and head 
vertex, then this is exactly the representation of the underlying undirected 
network. Thus, an efficient algorithm fcr finding cut-vertices in an undirected 


graph may be used to find cut-vertices in D. A slight modification to the 
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DFS (depth-first search) algorithm due to Tarjan [1972] can find all cut-vertices in 
O(|.£ |) time. 

Once all cut-vertices have been found and all irrelevant components 
removed, any cut-vertices remaining which are net K-vertices are neck vertices 
and may be added to the set K. It is fairly obvious that each component must work 
by itself if D is to work. Since edge reliabilities are independent it follows that 
Rex(D) = Rex(21) Rex (D2) ++: Pax(Dm) where the D; are the nonseparable, 
relevant components of J. This is identically equal to a network which has the 
same components as D but whose only cut-vertex is the source. Á raore rigorous 
demonstration of this identity can be made by applying the decomposition scheme 
described next, "moving edges to source," to all edges incident out of cut-vertices. 
Figure 3.2 demonstrates the results of cut-vertex decomposition applied to a 
network D with five components, two of which are irrelevant. In this case, 
Rex(D) = Rex(D1)*Rsx (Da) Fsx( Ds). 

2. Moving Edges to Source 

In this section, we discuss a technique which does not fit neatly under 
the rubric of decomposition or reduction since it only moves edges from one place 
to another in the network. However, by moving certain edge tails to the source s, 
new parallel reducticns, series reductions or contractions may be made possible 
along with cut-vertex decompositions. Consider an edge e=(u,u) such that vek. 
Edge e can be made incident out of any other K-vertex u'sK without afiecting 
reliability as long as no cycles are created. A simple conditioning argument can be 
used to show this, Agrawal [1982]. The source is the obvious candidate to receive 
the edge tails since it has no edges incident into it and therefore no cycles can 
possibly be created. Figure 3.3 shows an example where this operation creates a 


cut-vertex and allows a parallel reduction. 





IV. COMPUTING SAT XG O! 


This chapter gives a general method for computing SKT reliability in an acyclic 
directed network. The method is based on repeated computation of SAT reliability 
of certain subgraphs of D and, consequently, we first discuss SAT reliability. After 
deriving a method for computing SKT reliability, a simple lower bound becomes 
apparent and several heuristic arguments are given indicating that the bound’s 
accuracy should be good. Actual computational examples are reserved for the 


next chapter. 


A. SOURCE-TO-ALL TERMINAL RELIABILITY 


Let D(V,£) be an acyclic directed network and let g(2!—p, be the failure 
probability for any edge e;e£E and let E; be the set of edges incident into v,. Also, 
let the vertices be numbered such that (v,,v2, © : * Un} is an an acyclic ordering of 


V. Then, SAT reliability may be expressed as 


R,y(D) = P( s communicates with all veV—v, (4.1) 
and at least one edge into v, works ) 
since there no edges incident cut of v, by the acyclic ordering 


I 


> Rs(v-u Y (D—vn) 
9, بت‎ 


by independence 








A ADS 1] Ma H” 
rs 
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Another way to see this is to use the fundamental topological property that 
R,y(D) is unchanged if any edge (u,v) with veÄ is replaced by an edge (w,v) with 
the same reliability and weÄ has a lower number than u in an acyclic ordering of 
D. Consequently, in the SAT problem, we may move all edges so that they corae 
out of the source just like the "moving edge to source" operation in section 3.2. 
Equation (4.1) is then trivially true for a network where all edges come out of the 
source. (See Figure 4.1). The complexity of computing A,y(D) is easily seen to be 


O(|£ |) if all edges are with common head vertices are stored contiguously. 


B. SOURCE-TO-K-TERMINAL RELIABILITY 


Now we consider the SKT reliability problem in D. Let D(V,E) be a acyclic 


directed network with ACV. By simple state-space partitioning, 


Reg(D) = P{ s communicates with al vey) 


+ J P(s communicates with all veV—v, but not with v;) 
س ںا‎ kK 


+ >) P(s communicates with all veV—v,—v; but not with u; or uj) 
٭ رت پا‎ K 


E, .٦ 


+ P(s communicates with all veK but not with any ve V-K) 

But, 

P( s communicates with all veV) < R,y(D) 

P( s communicates with all v €V—v; but not with v ) 2 Fs(y-, (D). P( No edge into 
v; works ) 

P( s communicates with all v €V—v, —uj)7 Fs(y-4 (2). P( No edge into v; or v; 
works except an edge e z(vi,v;) or e z(v;,vi)) 

ت 


In general, then, 
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Rex(D) = Rey(D) + El (d, 0 2 = (4.2) 


اھ ور جا (A gR‏ + 
a‏ ۱ ]»» اج 
ug)‏ - وه - پره - )و2۱۳ ۲1 ]۲ + 
e) a ۱ ^ |‏ ور 
mw‏ 


where A; is the set of edges going into u,, Ay is the set of all edges going into v; or 
vj except edges going between v, and vj, Ej, is the set of all edges going into vi, vj 
or v, except edges going between pairs of vertices in (v.v; v, i, etc. 

The complexity of computing R,x(D) is, via the above formula, exponential in 


the number of non-X-vertices. The total number of additive terms in equation (4.2) 


n| 
is J CF = gI'-K! where C? is the number of ways to choose i elements out of n 
i=0 


elements. Computing SKT reliability thus requires O(|Z|2!V-*) time using 
equation (4.2). Of course, if the number of non-K-vertices is fixed as | F| changes. 


computation is effectively O(|£ |). 


C. A LOWER BOUND ON SKT RELIABILITY 


Since computation increases exponentially in the worst case, we may be 
satisfied with an approximation to SKT reliability which requires less time to 
compute than exact reliability. Since all terrns in equation (4.2) are positive, any 
subset of these terms provides an approximation to SKT reliability whica is a lower 
bound. Also, the product formed by lower bounds on independent components of a 
network yields a lower bound on overall network reliabiiity. The subset of terms we 
choose to form ¿his lower. bound consists of terms enumerated in the oraer given in 
equation (4.2) until some specified computational limit is reached. The motivation 


for this scheme being fairly accurate is given below. 
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Suppose p; = p for all e;eZ and p is close to 1. Consider the first few terms of 


Equation (4.2). 


12,1] 


R,y(D) = ah —Y سر‎ 0 


a du S.N 7 2 اقا ان‎ 
nu) q ani 1:8 q 


Ey] + | ۱و2 2 | (-| و5‎ 
. ۰ e < - re 

[E 85] ern? Ui U;) g md K q q 

Thus, the ccntribution of the terms to the total reliability sum diminishes 


rapidly unless there are a significant number of edges connecting non-K-vertices. 


Another reason that the lower bound should often be good is that terms with 
many vertices deleted may be identically zero. For example, consider an ST 
reliability problem where the shortest path from s to £ has l| edges in it. This path 
must contain | —1 non-K-vertices and so Fs (yu, =...) (Dv: —,—...) € 0 when 
|V-—hu Vj, . . . Ue S,6$|<l—1. Of course, other terms may be identically zero, too. 
For example, in any K-terminal problem, the term corresponding to Dv,—v, will 
be identically zero if all paths from s to some K-vertex v must include either v; or 
Uj. 

Comiputational experience in Chapter V lends credence to the accuracy of this 


lower bound. 





V. ALGORITHM AND DATA STRUCTURE 


The algorithm presented here for reliability evaluation of directed acyclic 
networks is based on the theory developed in chapter 3 and chapter 4. The 
program is intended for general use so it was coded in the widely available 


language, FORTRAN. This language is available on most computer systems and is 


still one of the most popular languages among operations research analysts. 


A | ALGORITHM 


The objective in this section is to develop an efficient algorithm for computing 
SKT reliability in a directed acyclic network. Initially, polynomial-time reductions 
and decomposition are used to reduce network parameters and decompose the the 
network into subnetworks, if possible. If reliability is not completely computed via 
reliability-preserving reductions, the general SKT reliability expression is 
computed as described in chapter 4. If a specified CPU time limit is reached for 
any component, the computation is truncated and a lower bound on reliability is 


produced. 


There is a one main algorithm and numerous subroutine algorithms in the - 
program. We discuss the main program and several important subroutine 
algorithms for understanding. The structure of the main program is based ona 
step by step approach which avoids redundant work. Comments will be given in 


curly brackets. 





Main 
Input: A directed acyclic network D with vertex set V, | V|z2, edge set E, | £22, 
and set KCV, | K|22. Edge reliability p; for each edge e; in £. 
Output: F,y(D), exact reliability, if. CPU time less than specified amount, otherwise 
a lower bound. 
1) {Initialize} M = 1.0, bound=.false. 
2) Read all edge data and create network data structure 


3) {Change neck vertices to K-vertices} K = Neck (D,K) 


( 
( 
( 
(4) {Perform reductions for whole network? Reduct (D,K,M) ° 
(5) If remaining vertices | V| < 1, go to end. 

(6) §D,,Do,....D{= Decomp(D) 

(7) For each subnetwork J; 

(a) {Initialize} lowbnd=.false., set time = 0. 

(b) {Move edge tails to source} For each edge e=(v,,u;) with ue 
{move v, tos} Hove(e) 

(c) Construct the data structure for subnetwork D, 

(d) {Reduce subnetwork] Reduct (D, K;, M) 

(e) § Initialize } sumpro=0 

(f) { Find SAT reliability } sumpro = Rsa(D,), 

(g) Count the remaining non-K-vertices for D,, i.e., if vertex number 


> 0, count=count+i 


1) {Find the next combination of non-K-vertices} 
tuj, Uj,.. uU] 7 Comb( i) 
2) Í Find ilgi for رت‎ di ! M' = Multi(v; .;,. NS L7» 


3) sumpro=sumpro + MH" Rsa(D—v; زیت‎ =...) 


2 


4) If elapsed CPU time = time limit then lowbnd=.true. and 
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bound-.true. 


5) If lowbnd=.false. and all combination have not been 
enumerated then go to (B.h.1) 
(i) M=M x sumpro 
(8) If bound=.false. print "Exact reliability is" 4, otherwise print "Lower 
bound on reliability is" M 


End Main 


function Neck (D,K) 
Input: D,K 
Output: New set of It-vertices 
( This routine finds all neck vertices and changes them to K-vertices ) 
(1) Put all v, K into Q, 
(2) cz | K | 
(3) While Q, £2 
(a) Remove v from Q, 
(b) Search all K-vertices in D—v, using Depth-l'irst Search 
(c) k'z number of K-vertices reached 
(c) If &'«k then K 2 KUJfu] 
(3) Return (A). 
End of Neck 


subroutine Reduct(D,K,M) 

Input: D,K,M 

Output: Reduced network D with modified K and /4 
( This routine performs all non-K-vertex reductions, series-parallel 


reductions and indegree-one K-vertex contractions ) 


(1) Construct Q; 2$ veV-K | outdeg(v)-0 | 
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(2) While Q, ~¢ 
(a) For each vertex v; incident to wu, if outdeg(v;)=1 and vjeV—-K 
then add v; to Q4. 
(b) Nonkre ( 1) 
(3) Construct Q, - | weV-K | indeg(v;¡)= outdeg(v;)-1 ] and & = | wei | 
indeg(vi)71 | 
(4) For each v, 
(a) Check all outgoing edges, if e(u,vj),e(w,v;) are found then 
Parall( e) 
(b) If parallel reductions creates new series or contraction vertices, 
put these into Q, or Q., respectively. 
(5) While Q,#¢ and Q, #¢ 
(a) lf Q.4#¢ Series( Q, (1)) 
1) If series reductions create new paraliel edges, do parallel 
reduction. Parall( e) 
2) If parallel reductions create new series or contraction 
vertices, put these into Q, or Q., respectively. 
(b) If Q; 4 Contra( Q. (1)) 
1) If the contraction creates a new contraction vertex, put 
these into Q. 
2) If the contraction creates new parallel edges, do parallel 
reduction Parall(2 ) 
3) If parallel reductions create new series or contraction 
vertices, put these into 9, or @,, respectively 


(6) Return (D,K,M) 


End Reduct 





subroutine Decomp(D) 
Input: D 
Output: Subnetworks D,,D2...., Dm 
(1) While not all cut-vertices have been found 
(a) Using DFS, find pendant component D, with cut-vertex v, 
(b) If Z2; has no K-vertices except possibly v, then delete D; 
(9) Return (2, Ds..... Eua) 
End Decomp 


After reading in the data and creating the network data structure, the 
algorithm begins by changing any neck vertices to K-vertices by a call to Neck at 
step (3). Next, a call to Reduct is made at step (4). Reduct first removes all 
outdegree zero non-K-vertices. Through the queue mechanism, it is ensured that 
all initial outdegree zero non-K-vertices are deleted along with any which are 
created as the reduction proceeds. Á similar method is used for the other 


reductions schemes. 


The other reductions, indegree-one-K-vertex contractions and parallel and 
series reductions, may recursively create new reducible vertices or edges. The 
newly created reducible vertices are put into a queue instead of immediately being 
reduced to avoid redundant work. Parallel edges are reduced immediately, 
however. The network of Figure 5.1 is completely reduced by reductions alone and 
illustrates how a reducticn can create other reducible edges or vertices. A call to 
Reduct would result in the set of reductions given below. This network would be 


completely reduced and its reliability completely calculated by these reductions. 





After 


(1) Reduce parallel edges 
(2) contract vertex 
(3) series vertex 
(4) contract vertex 
(5) reduce parallel edges 
(6) series vertex 

(6) series vertex 

(7) reduce parallel edges 
B) series vertex 

9) reduce parallel edges 


10) contract vertex 


12) contract vertex 


( 
( 
( 
(11) reduce parallel edges 
( 
(13) reduce parallel edges 
( 


14) contract vertex 


reductions are 


performed, a 


6۵,6 0 و‎ 
Ug 

Va 

Uz 

85,85 LO و‎ 
v? 

Ur 

6 ۱0,6 ۱۱ ۲۵ € 1g 
Ug 

€ 10.8 12 LO € 1g 
Ug 

92,6 و‎ LO 85 
Ug 

84,9045 LO 23 


Us 


network will 


subeomponents if it contains any cut-vertices. This is done at step (6) of Main bya 
call to Decomp. The rest of the calculations are performed iteratively for each 
component D, under step (7). First, all edges with tail-vertices in X are moved to 
the source. This may produce additional reductions, so another call to Reduct is 
made for the subcomponent. Then, the terms corresponding to Equation 4.2 (for 
that component) are produced until all such terms have been enumerated or the 
CPU time limit is reached. Note that this CPU limit is based on the CPU time used 


in calculating subcomponent reliability and is not based on total CPU time. 


If the CPU time limit is reached for any of the subcomponents, the final value 
obtained by the algorithm is a lower bound on reliability. Otherwise, the value 


obtained is exact reliability. Choosing a CPU time limit will depend on computer 
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facilities available and the network under analysis. Figure 6.5 shows the reliability 
increment as a function of CPU time. The lower bound obtained using a reasonable 
CPU time limit will not differ much from the exact reliability for this large scale 


network example 


B. DATA STRUCTURE 


There a many ways to represent a network on a digital computer. The most 
straightforward is to use an adjacency matrix: If. vertices v; and v; are adjacent, 
M(i,j)=1, and M(ij)-0 otherwise. Matrices of the above type are, in practice, 
extremely sparse and this representation is very inefficient in terms of space. ít is 
often ineflicient in terms of execution time, too, particularly when one is 
interested in operations that require retrieving information from all vertices 
adjacent to a particular vertex, e.g., Neck, Decomp, Reduct, and Rsa. In this 
situation, since most of the elements of M(ij) are null, the algorithm will spend a 
great amount of time retrieving and comparing zero values. For a static network, 
an efficient representation is a packed matrix. In our algorithm, the network is 
dynamic, since it is reduced and decomposed many times during execution. 
Therefore, a multi-linked data structure (multiple linked lists with links to other 
data structures such as arrays) is used to represent the network. Every vertex has 
a linked list of adjacent vertices, which, besides informing which vertices are 
adjacent to it, also tells whether or not the vertex belongs to set K. Two such lists 
are kept for each vertex, an “adjacent out” list and an "adjacent into” list. OF 
course, each adjacency entry corresponds to an edge. So, in addition, there are 


pointers indicating the addresses where information about edges is Kept. 


Figure 5.3 illustrates this data structure. The vertex v, only can send 


communication to vertices vo, and va, and vertex va can receive the 





communication from vertex v, and send the communication to vertices v4 and v,, 


and so on. This structure is a very efficient data structure for dynamic networks. 


C. ALGORITHMIC COHPLEXITY 


It is important to determine the complexity of the various parts cf the 
algorithm so that overall complexity can be understood. In this section, we analyze 
the worst-case complexity of the subroutines, as written, and thus, find the 
complexity of the overall algorithm. Not all of the subroutines have been written 
as efficiently as possible because cf the difficulty in programming such routines 
and because of limited time. Furthermore, for problems of moderate size, it is 
unlikely that much efficiency would be gained using more sophisticated routines, 
since, in practice, most of the routines written seem to operate mere efficiently 


than the worst-case analysis indicates. 


There are three different reduction schemes included in the algorithm: non- 
K-vertex deletion, indegree-one K-vertex contraction and series and parallel 


reductions. However, the last three of these interact directly. 


The non-K-vertex deletion, steps (1) and (2) in REDUCT, requires O(|£Z |) 
operations in the worst case. For example, consider a complete acyclic graph 
which has only one K-vertex, the source. Here, one vertex at a time would oe 
deleted from the network, but, in the process, every edge would be examined 
exactly once until a single isolated vertex remained. The initial queue-building is cf 


no consequence since it is O(| Vl). 


Steps (3) through (5) constitute the other three reductions. Queue-building 
again may be disregarded since it is at an O(!V|)} operation. Initial parallel 
reductions may require O(| Z |) operations using a bucket sort technique to check 


for edges with common end vertices. The central step of the algorithm, step (5), 





requires O(|V|*) time. This is true since (1) at most |V| series reductions or 
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contractions can ever occur (2) each such reduction requires one O(| V|) check for 
a newly created parallel edge, and (3) once identified, actually carrying out the 
parallel reduction requires only constant time. Overall then, REDUCT requires 
O(| V|*) time. 

There are two different decomposition schemes used, cut-vertex 
decomposition and neck-vertex decomposition. The complexity of the cut-vertex 
algorithm is O([E|]) since it is based on the O(|Z|) depth-first search. The 
complexity of the neck-vertex algorithm will be O(|V—-K||E |) since finding neck- 
vertices is based on repeated O( | Z |) searches on D—v for each non-K-vertex v. in 
the worst case then, this decomposition requires O(| V! | Z|) time. However, if we 
fix the number of non-K-vertices as described in chapter 4, then the complexity of 


neck-vertex algorithm will effectively be O( | Z |). 


It follows from the above discussion that tbe total complexity of the reduction 
and decomposition algorithms is O(|V|1Z |) nermally, and O(|V|*) when non-K- 


vertices are fixed. 


After reductions, computing SKT reliability requires O(! £ |2!V-Kl) time. This is 
obviously true since 2!V-Xl combinations of non-K-vertices must be produced and 
an O(| Z |) SAT computation made for each combination. When the number of non- 
K-vertices is fixed, the computation becomes O(|£'|). The entire algorithm 
therefore requires O(|£ |21Y-X1) time in general but only O(|V!?) time with non-K- 


vertices fixed. 
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VI. COMPUTATIONAL EXPERIENCE, 


In order to test the efficiency of the SKT reliability algorithm, we used two 
kinds of networks, “complete” acyclic networks, and acyclic "street" networks. A 
complete acyclic network is a graph D(V,£) such that V=fu,,V2,*** ,Unj and 


e-(w,vj)eE if and only if i<j. Any such network with n vertices has exactly 


۳۳ edges A network is called an m xn "street" network if (1) its vertices can 


be laid out in an m rows by n columns rectanguiar grid, (2) for all vertices except 
those in the last column there is an edge incident out each vertex directed to the 
right, and (3) for all vertices except those in the first row, there is an edge incident 
out of each vertex directed upward. When m 271 case the street retwork is calied 


"square" (See Figure 6.1). 


In the complete network we assigned v, to be the source and fixed v, as a K- 
vertex. This insures that all edges and vertices are relevant. Other K-vertices were 
chosen randomly as required. In the street network, all corner vertices were fixed 
as K-vertices, the lower left-hand corner vertex assigned as the source and all 


other K-vertices chosen randomly. 


A | SOURCE-TO-ALL-TERMINAL RELIABILITY COMPUTATION 


As described in chapter 4, the complexity of computing SAT reliability is 
O(|E]|). Since repeated computation of SAT reliability is the key to computing 
generai SKT reliability, it is interesting to see how much CPU time is required for a 
single SAT computation. ligure 8.2 shows the CPU time used as a function of the 
number of vertices in the two types of test networks. |V| is increased from 15 to 
80 for the complete network and from 9 to 144 for the square street network. Of 


course, es a function of !V!, | Z! is quadratic for the complete network ranging 
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VL  COMPUTATIO 


In order to test the efficiency of the SKT reliability algorithm, we used two 
kinds of networks, "complete" acyclic networks, and acyclic "street" networks. A 
completa acyclic network is a graph D(V,E) such that V-ívu,,vs,-: - Un} and 


e-(w,vj)eE if and only if i<j. Any such network with n vertices has exactly 


n(n-1) n= edges. A network is called an m xn "street" network if (1) its vertices can 


be laid out in an m rows by n columns rectanguiar grid, (2) for all vertices except 
those in the last column there is an edge incident out each vertex directed to the 
right, and (3) for all vertices except those in the first row, there is an edge incident 
out of each vertex directed upward. When n 27 case the street network is calied 


"square" (See Figure 6.1). 


In the complete network we assigned v, to be the source and fixed v, as a K- 
vertex. This insures tbat all edges and vertices are relevant. Other K-vertices were | 
chosen randomly as required. In the street network, all corner vertices were fixed 
as K-vertices, the lower left-hand corner vertex assigned as the source and all 


other K-vertices chosen randomly. 


A. SOURCE-TO-ALL-TERMINAL RELIABILITY COMPUTATION 


As described in chapter 4, the complexity of computing SAT reliability is 
O(|£|). Since repeated computation of SAT reliability is the key to computing 
generai SKT reliability, it is interesting to see how much CPU time is required fora 
single SAT computation. Figure 6.2 shows the CPU time used as a function of the 
number of vertices in the two types of test networks. |V| is increased from 15 to 
80 for the complete network and from 9 to 144 tor the square street network. Di 


course, as a function of |V!, | E! is quadratic for the complete network ranging 
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from 105 to 3160. |Z| increases linearly as a function of |V| for the street 
network ranging from 12 to 2B4. From the figure, it is clear that SAT reliability can 


be cemputed very rapidly. 


B. SOURCE-TO-K-TERMINAL COMPUTATION 


Here we investigate the work required by the algorithm for computing general 
SKT reliability and the efficiency gained by using reductions and decomposition. 
Initial reductions and decomposition require O(|V| | £]) time. After reductions and 
decomposition, the complexity of | computing OKT reliability is 


O(maxt (A ||£,|3) where C is the set of separable components of D and where 
[13 


VW -K; is the set of non-K-vertices in component i. Consequently, we expect an 
exponential growth in the worst case for SKT computations. Figure 6.3 shows tne 
growth in CPU time as the number of non-K-vertices is increased, with and withcut 


reductions and decomposition. 


For the complete network, the number of non-K-vertices ranged from 3 to 14 
out of a total 20 vertices. The reductions and decomposition do not reduce 
computation much here since the complete network is so densely connected. Un 
the other hand, computation times are significantly reduced when reductions and 
decomposition are used in the street network, primarily because many non-K- 
vertices are changed to K-vertices. Notice that the computation time required bv 
the street network without reductions and decomposition is almost the same as for 
the complete network. This shows how complexity is dependent mostly on tne 


number of non-K-vertices. 


C. COMPUTATION WITH NON-K-VERTICES FIXED 


If the number of non-K-vertices is fixed, ie. |V-A|=c for some constant c, 
then the complexity of computing SKT reliability is O(| V|)-0(c | E]) - O(2* E) 


which is O(|V|?). Thus, the complexity is polynomial instead of exponential. Of 
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course, if reductions and decomposition are not used, the complexity becomes 
00۱۳ ۱( which is better than O(|V|*). However, in practice, the reductions and 
decomposition are very efficient and do not seem to increase computation time. 
Figure 6.4 shows the increase in computation for both test networks as |V| is 


increased. 


D. LOWER BOUND COMPUTATION 


As described in cnapters 4 and 5, a lower bound on SKT reliability is obtained 
when the normal SKT reliability computation is truncated. Two heuristic 
arguments were given in chapter 4 indicating that this tower bound should be fairly 
good. Figure 6.5 shows, for actual computations, how the lower bound approaches 
the exact network reliability very quickly. The test networks were assigned 
different values of p, to show that the accuracy of tne bound is not much affected 
by nı. It appears that a CPU time limit can be used successfully for truncating SKT 


computation to obtain a good lower bound when exact computation is not feasible. 





VII. CONCLUSION 


This thesis has explored the source-to-K-terminal reliability analysis problem 
in acyclic directed networks. We have shown that the polynomial-time reductions 
and decompositions reduce an exponentially complex problem to a polynomially 
complex problem under certain conditions. Even when these conditions do not 
hold, surprisingly large networks can be analyzed. The theoretical and 
computational results obtained here can be used by researchers in analyzing the 
reliability of directed acyclic networks and, in some cases, for the partial 


reliability analysis of general directed networks. 


We assumed that all vertices work perfectly and that edge reliabilities are 
constant while the network system is alive. These assumptions may not be 
acceptable in a real-world situation. It may be necessary to treat vertex 
reliabilities explicitly and. in some way, handle component reliabilities changing as 
a function of time. This will complicate reliability analysis but does not make 


analysis infeasible. 


The extension to unreliable vertices is nol hard. All K-vertices must work if 
the system is to work and thus, for unreliable K-vertices, we obtain network 
reliability as usual and then multiply by the probability that all K-vertices function. 
Unreliable non-K-vertices may be handled in several ways. One of the simplest 
ways is to split each unreliable non-K-vertex into two perfectly reliable vertices 
with an unreliable edge between the two vertices with a failure probability equal to 
the failure probability of the original vertex. All edges incident into the original 
vertex are made incident into the tail vertex of the new edge and all edges incident 


out of the original vertex are made incident out of the new head vertex. However, 
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it is also possible to redefine the reduction and decomposition schemes and 
redefine tbe general reliability expression to directly handle unreliable non-K- 
vertices. For example, the series reduction can include the ver*ex reliability so 
that PD, = pape p (",). Also, the terrm in the general reliability expression associated 
with "no edges in Ay. work" can be replaced with.a term wnich expresses "no 
edges into v, work or v; fails, and no edges into v; work or v; fails, and ... except 
edges going between pairs of vertices ih (9p7; cuj] This is not hard we 


compute. 


Reliability changing as a function of time can be hendled in one of two ways. 
The easiest way is to iust input component reliabilities at different times and 
compute network reliability at these times. Interpolation can be used to 
approximate values between the times for which reliability is actuaily evaluated. 
Another possibility is to produce a symbolic reliability expression that can be 
repeatedly evaluated by a computer at different times, thus avoiding repeated 
reductions, decompositions, etc. This might require large amounts of storage but, 
if reliability can be calculated at all in a reasonable amount of time, it should be 
feasible. Once the price is paid, computing reliability for any given point in time 


should be very eflicient. 


As modern computer, power and other network systems become larger and 
more complicated, reliability becomes a more important part of design criteria. 
Unfortunately, it is obvious that as networks become more complicated they also 
become more difficult to analyze. Further research should concentrate on ways to 
streamline computations and to find accurate but easily computable bounds. The 
method described in chapter 4 for computing SKT reliability could be improved if 
only non-zero terms were generated; An enumeration scheme similar to that of 
Satyanaryana and Prabhakar [1978] should help in this respect. The lower bound 


described in this thesis is a good start on finding accurate and easily computable 





bounds, but every lower bound shouid be accompanied with an upper bound so that 
accuracy can be checked. Finding a good upper bound should be a priority for 


future research. 
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APPENDIX A 
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V à Cut-vertex 


Bar, 





u 
j (b) ۰ 
v not a cut-vertex 
Neck Vertex to K-vertex Reduction 
«Figure 3.1> 
Decomposition 
<Figure 3.2> 
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Moving Edges to Source 
«Figure 3.3» 


h 
| 5 
y 
2 E 
S : oa, ==> s B 
E Qv 
3 


2 y 


All Terminal Reliability ",v(Di) = Rey (De) 
«Figure 4.1» 





V 


e, 13 2 14 o? 
6 
~é 
OL ss 
رت کو سد‎ 
| 7 





وھ ر 


سے یگ 


V4 


Comp!etely Reducible Network by Reduction Schemes 
«Figure 5.1» 





Ad). out 


Multilist Data Structure 
<Fig ace 9.2> 
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4x4 Square "Street" Network 
«Figure 6.1» 





Source-To-All Terminal Network 
< Figure 6.2 > 
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ELAPSED TIME: SEC. 
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< Figure 6.3 > 


42 
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REILABILITY INCREASE BY TIME IN TOTAL 20 VERTICES 
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< Figure 6.5 > 
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